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Ice  Forces  on  Flat,  Vertical  Indentors 
Pushed  Through  Floating  Ice  Sheets 

NAOKI  NAKAZAWA  AND  DEVINDER  S.  SODHI 


INTRODUCTION 

About  half  the  Earth’s  surface  is  subject  to  snow,  ice 
and  seasonally  frozen  ground.  Extensive  permafrost 
and  ice  sheets  exist  in  the  polar  regions,  and  the  oceans 
are  covered  by  seasonal  and  multiyear  sea  ice.  With  the 
growth  of  human  activities  in  cold  regions,  the  impact 
of  snow,  ice  and  frozen  ground  on  human  affairs  and 
commercial  development  has  grown  in  recent  decades, 
thus  establishing  the  need  for  scientific  research.  In  low 
temperatures,  the  properties  of  most  materials  change, 
and  many  machines  either  do  not  work  as  designed  or 
fail.  Design  criteria  and  construction  techniques  for 
structures  for  cold  regions  are  different  from  those  in 
warm  regions. 

Purpose  for  this  study 

A  rational  basis  is  needed  for  the  design  and  con¬ 
struction  of  structures  in  coastal  and  offshore  regions 
where  the  presence  of  floating  ice  presents  a  hazard. 
Among  the  ice  forces  imposed  on  structures  that  interact 
with  ice,  the  impact  of  an  ice  floe  produces  some  of  the 
greatest  loads  that  the  structure  has  to  be  designed  to 
withstand.  This  impact  can  be  approximated  by  an 
indentation. 

In  this  study,  indentation  tests  were  conducted  to 
allow  us  to  observe  ice  failure  and  to  determine  the 
distribution  of  ice  pressures  and  the  total  ice  force.  The 
geometry  of  an  indentation  in  an  actual  environment  is 
shown  in  Figure  1 .  An  ice  sheet  of  thickness  It  moves  at 


a  velocity  v  past  an  indentor  of  width  </.  Because  of  the 
limits  imposed  by  model  testing  in  the  laboratory, 
indentors  with  a  velocity  v  were  pushed  through  an  ice 
sheet  in  this  study. 

Previous  indentation  studies  were  conducted  using 
finite-sized  ice  sheets  that  were  confined  in  a  frame. 
Because  of  the  limitations  of  hydraulic  equipment  or 
facilities,  many  of  these  tests  used  a  short  indentation 
distance.  However,  limited  studies  by  Kato  and  Sodhi 
(1984),  Sodhi  and  Morris  (1984)  and  Tirnco  (1987) 
have  been  conducted  on  large,  floating  ice  sheets  with 
long  indentation  distances,  simulating  the  indentation 
of  a  vertical  structure  intoan  infinite,  first-year,  floating 
ice  sheet. 

The  previous  studies  did  not  completely  illuminate 
the  ice  failure  mechanism.  It  is,  therefore,  the  objective 
of  this  study  to  attempt  to  understand  the  ice  failure 
process  during  ice  crushing  against  an  indentor.  To 
accomplish  this,  instrumented  indentors  were  pushed 
against  the  edge  of  floating  ice  sheets  at  different 
velocities.  The  instrumentation  of  the  indentors  al¬ 
lowed  us  to  measure  both  the  total  ice  force  and  the  local 
ice  pressure  or  the  position  of  the  resultant  ice  force 
within  the  contact  area.  Furthermore,  transducers  were 
placed  to  monitor  the  Acoustic  Emission  (AE)  activity 
in  the  ice  and  to  measure  displacements  of  the  carriage 
and  the  indentor,  which  enabled  us  to  analyze  the 
energies  stored  in  the  structure  and  dissipated  in  the  ice. 

Plots  of  force  versus  time  and  force  versus  displace¬ 
ment  allowed  us  to  draw  conclusions  concerning  the 


Figure  1.  Schematic  of  indentation  geometry. 


magnitude  of  the  force  or  energy  required  to  make  an  ice 
sheet  fail,  as  well  as  allowing  us  to  observe  the  size  of 
the  damage  zone  in  the  ice  sheet  during  each  loading 
event.  Moreover,  AE  signals,  believed  to  be  caused  by 
the  formation  of  microcracks,  were  correlated  with  the 
damage  caused  by  the  interaction  and  with  the  resulting 
ice  forces.  To  summarize  the  objectives  of  this  study, 
we  hoped  to  do  the  following: 

1 .  Observe  the  nature  of  the  ice  failure  mode  during 
the  indentation. 

2.  Estimate  the  energy  required  to  cause  an  ice  sheet 
to  fail. 

3.  Observe  the  indentor  response  at  the  point  of  the 
ice  failure. 

4.  Determine  whetherthere  is  any  non-simultaneous 
crushing  failure  by  measuring  the  distribution  of  ice 
pressure  at  the  ice/structure  interface. 

5.  Determine  the  frequency  of  ice  crushing  failure. 

The  instrumentation  used  in  this  study  had  some 

desirable  features.  First  of  all.  the  experiments  were 
conducted  on  large,  floating  ice  sheets,  simulating  an 
infinite  ice  sheet;  second,  the  indentorsupport  wasquite 
stiff;  third,  a  screw-driven  carriage  capable  of  moving 
up  to  2  m  was  used;  fourth,  direct  ice  force  measure¬ 
ments  were  made  by  mounting  the  indentor  on  the  load 
cells  at  the  ice/structure  interface;  and  fifth,  separate 
measurements  were  made  of  carriage  and  indentor 
displacement  relative  to  an  ice  sheet. 

We  carried  out  92  indentation  tests  with  rigid,  verti¬ 
cal,  flat  indentors  of  various  widths  at  different  inden¬ 
tation  velocities  on  different  thicknesses  of  freshwater 
ice.  We  varied  several  parameters  during  this  study — 
indentor  widths  of  50. 60,  100  and  150  mm  were  used, 
ice  thickness  ranged  from  20  to  60  mm,  and  indentor 
velocity  varied  between  1  and  9  mm/s. 

During  each  test,  we  measured  the  force  generated 


during  ice-indentor  interaction,  ice  pressure  at  the  ice/ 
indentor  interface,  indentordisplacementrelativetothe 
ice  sheet,  displacement  of  the  carriage  relative  to  a  fixed 
datum,  and  AE  activity  in  the  ice.  In  addition,  we  noted 
nucleation  of  the  first  radial  crack  using  an  event  marker, 
and  used  photography  and  video  to  record  the  experi¬ 
ments. 

Background 

The  crystallographic  structure  and  properties  of  both 
freshwater  and  sea  ice  have  been  intensively  studied 
(e.g..  Weeks  and  Ackley  1982,  Weeks  and  Cox  1984). 
Through  studies  by  many  investigators,  the  dependence 
of  ice  strength  (i.e.,  compressive,  tensile  and  bending 
strength)  on  temperature,  strain  rate,  grain  size,  poros¬ 
ity,  salinity,  etc.,  is  now  well  understood.  Frictional 
resistance  and  adfreeze  bond  strength  between  ice  and 
different  materials  (e.g.,  concrete,  steel  and  wood)  have 
been  investigated  by  Tusima  and  Tabata  ( 1 979),  Oksanen 
( 1 980),  Forland  and  Tatinclaux  ( 1 984)  and  Saeki  et  al. 
(1986,  1988).  Ice  forces  on  bridge  piers  in  rivers  have 
been  studied  for  a  long  time,  but  the  interaction  between 
ice  and  structures  in  offshore  regions  is  a  relatively  new 
field  of  study. 

For  the  estimation  of  design  loads,  ice  forces  are 
broadly  divided  into  two  categories  (Neill  1976,  San¬ 
derson  1988):  static  and  dynamic. 

Static  loading 

The  loading  state  is  defined  as  static  if  ice  exists  in 
stationary  contact  with  a  structure,  and  then  the  struc¬ 
ture  experiences  an  increasing  load  applied  to  it  by 
natural  driving  forces,  such  as  wind  and  water  stresses, 
and  thermal  expansion  of  ice  sheets  because  of  warm¬ 
ing. 

Dynamic  loading 

The  loading  state  is  defined  as  dynamic  if  an  ice 
feature  is  not  initially  in  contact  with  a  structure,  but 
arrives  and  strikes  it  with  an  appreciable  velocity.  This 
loading  state  differs  significantly  from  static  loading  in 
two  respects:  firstly,  the  initial  contact  conditions  are 
invariably  irregular  and  nonuniform,  and  secondly,  the 
duration  of  the  impact  is  generally  determined  by  the 
kinetic  energy  of  the  impacting  ice  feature,  which  may 
come  to  rest  during  the  process.  Examples  of  dynamic 
ice  forces  are  impacts  by  multiyear  floes  and  by  ice¬ 
bergs  or  ice  islands. 

In  addition  to  the  above  two  broad  categories  of 
loading,  it  is  also  necessary  to  understand  the  magnitude 
of  the  forces  and  the  contact  area  over  which  these 
forces  act.  For  example,  total  force  on  the  whole  struc¬ 
ture  and  local  pressures  over  a  limited  area  are  ex¬ 
tremely  important. 


Total  load 

The  total  load  sustained  by  a  structure  is  important 
for  considerations  of  foundation  sliding  resistance, 
foundation  bearing  capacity  and  overturning  moment. 

Local  loads  or  pressure 

The  magnitude  and  distribution  of  local  pressure  is 
essential  information  for  determination  of  the  design 
and  spacing  of  internal  structural  members  and  the 
dimensions  of  internal  cell  units.  Loads  over  smaller 
areas  are  also  important  for  the  design  of  the  external 
skin  of  a  structure. 

Structures  placed  in  an  ice  environment  should  be 
able  to  withstand  not  only  the  total  ice  load  but  also  the 
local  pressure.  A  good  understanding  of  these  loads  will 
lead  to  an  economical  design  for  structures. 

Experimental  and  theoretical  studies  on  ice  pressure 
distribution  at  the  ice/structure  interface  were  con¬ 
ducted  by  Schwarz  ( 1 970),  Kry  ( 1 978, 1 979),  Tanaka et 
al.  ( 1 987),  Tunik  (1987)  and  Blanchet  ( 1 987). 

Korzhavin’s  formula,  which  empirically  relates  ice 
crushing  force  with  contact  area  and  compressive 
strength  of  ice,  has  been  discussed  by  many  researchers. 
Experimental  and  theoretical  studies  to  determine  ice 
forces  on  structures  were  conducted  by  many  research¬ 
ers:  laboratory  tests  were  done  by  Hirayama  et  al. 
( 1 974),  Michel  and  Toussaint  ( 1 977),  Saeki  et  al.  ( 1 977), 
Kry  (1980b),  Michel  and  Blanchet  (1983),  Kato  and 


Sodhi  (1984),  Sodhi  and  Morris  (1984)  and  Timco 
(1987);  field  tests  have  been  done  by  Zabilansky  et  al. 
(1975)  and  Croasdale  et  al.  (1977). 


EXPERIMENTAL  SETUP  AND  PROCEDURES 
Facilities 

The  experiments  were  conducted  in  the  test  basin  of 
CRREL's  Ice  Engineering  Facility.  The  test  basin  is 
34.4  m  (1 13  ft)  long.  9  m  (30  ft)  wide  and  2.4  m  (8  ft) 
deep.  A  photograph  of  the  test  basin  is  shown  in  Figure 
2.  The  test  basin  is  insulated,  and  its  refrigeration  is 
provided  by  seven  forced-air  heat  exchangers  suspended 
from  the  ceiling.  The  total  refrigeration  capacity  is  80 
kW,  and  the  minimum  air  temperature  that  can  lie 
achieved  is  approximately  -23°C.  The  refrigeration 
fluid  is  ammonia. 

Ice  sheets 

Ice  growth 

The  freshwater  ice  sheets  were  grown  in  the  test 
basin.  To  freeze  an  ice  sheet,  a  mixture  of  water  and  air 
was  sprayed  i nto  the  cold  air  over  the  water  surface.  The 
resulting  seed  crystals  falling  from  the  air  onto  the  water 
surface  started  the  ice  growth  and  ensured  uniform, 
small  size  grains  (approximately  1  to  2  mm)  throughout 
the  sheet.  The  top  layer  of  the  ice  sheet  is  generally  thin 
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( 1  mm)  and  is  composed  of  columnar  ice  with  a  mixed 
vertical  and  horizontal  c-axis  orientation.  The  rest  of  the 
ice  sheet  has  a  random  c-axis  orientation  in  the  horizon¬ 
tal  plane  and  has  a  columnar  structure  characteristic  of 
dendritic  type  growth. 

The  ice  sheets  were  grown  at  several  temperatures, 
depending  on  the  time  available  for  freezing  and  the 
cooling  capacity  available  from  the  refrigeration  plant. 
Most  sheets  were  grown  at  a  temperature  between  -1 5 
and-20°C.  It  took  approximately  36  hours  to  grow  a  50- 
mm-thick  ice  sheet.  On  a  typical  day  of  testing,  the  room 
was  allowed  to  wami  to  approximately  -3°C,  after 
which  the  ice  was  allowed  to  temper  and  to  attain  a 
uniform  temperature. 

Temperature 

The  temperature  of  the  ice  sheets  was  measured  a 
few  times  during  the  program  and  was  found  to  be 
between -1  andO°C. 

Ice  thickness 

The  ice  thickness  was  measured  after  each  test  at  the 
test  track  using  vernier  calipers  with  a  resolution  of  0.05 
mm. 

Characteristic  length 

The  characteristic  length  of  each  floating  ice  sheet 
was  measured  shortly  before  the  tests  for  that  day.  A 
deadweight  (10  lb  [5  kg])  was  placed  over  a  circular 
plate  near  the  center  of  an  ice  sheet.  The  resulting 
deflection  of  the  sheet  was  measured  using  a  displace¬ 
ment  transducer  at  the  point  of  loading.  The  character¬ 
istic  length  of  the  floating  ice  sheet  was  calculated 
according  to  the  theory  of  an  infinite  plate  on  an  elastic 
foundation  (Wyman  1950,  Sodhi  et  al.  1982). 

The  setup  for  load  application  and  deformation 
measurements  is  shown  schematically  in  Figure  3.  The 
characteristic  length  was  calculated  using  eq  1 . 


Ice  Sheet 


Figure  3 .  Experimental  setup  for  measuring  the  charac¬ 
teristic  length  of  a  floating  ice  sheet. 


where  /  =  characteristic  length  of  ice  sheet 

A P  =  incremental  load  placed  on  the  ice  sheet 
K  =  specific  weight  of  water 
Aw  =  resulting  deflection  increment  of  the  ice 
sheet  at  the  center  of  the  load  zone 
lny  =  0.55772157  (Euler  constant) 
r  =  radius  of  the  applied  load 
a  =  ill. 

Elastic  modulus 

The  elastic  modulus  can  then  be  calculated  by  eq  2 
using  the  characteristic  length  and  ice  thickness 

r  \2{\-v2)K  /4 

t  =  -  iii 


where  £  =  effective  elastic  modulus 
v  =  Poisson's  ratio  for  ice 
K  =  specific  weight  of  water 
/  =  characteristic  length  of  the  ice 
h  =  ice  thickness. 

Experimental  apparatus 

A  schematic  drawing  of  the  test  structure  setup  is 
shown  in  Figure  4.  A  carriage  that  travels  parallel  to  the 
length  of  the  basin  was  used  to  carry  the  test  apparatus. 
The  high-force  module,  which  includes  a  screw-driven 
carriage  for  the  indentation  tests,  was  mounted  under 
the  carriage.  The  test  structure,  structural  support  and 
indentor  were  attached  to  the  screw-driven  carriage. 
The  test  structure  traveled  perpendicular  to  the  longitu¬ 
dinal  direction  of  the  basin,  and  was  powered  by  a  motor 
located  at  one  end  of  the  high-force  module.  The  range 
of  velocity  for  the  carriage  in  the  high-force  module  was 
from  1  to  9  mm/s.  A  photograph  of  the  high-force 
module  is  shown  in  Figure  5. 

A  schematic  diagram  of  the  experimental  setup  is 


Control 

Room 


Figure  4.  Schematic  of  the  test  structure  setup. 


shown  in  Figure  6.  The  structural  support  mounted  on 
the  carriage  of  the  high-force  module  is  1 .2  m  long.  0.6 
m  wide  and  0.3  m  high.  The  indentor  support  was 
mounted  under  the  structural  support,  and  indentor 
plates  were  attached  to  the  indentor  support.  Both  the 
indentor  and  structural  support  moved  together.  The 
load  cell  installed  in  the  structural  support  was  cali¬ 
brated  to  measure  the  horizontal  ice  force  on  the  inden¬ 
tor. 

The  following  two  methods  were  used  to  determine 
the  pressure  distribution  on  the  indentor  (Fig.  7). 

Installation  of  pressure  transducers 
on  the  indentor  plate 

The  pressure  transducers  were  installed  to  measure 
the  local  ice  pressure  at  the  ice-indentor  interface.  The 
widths  of  the  indentor  used  with  this  method  were  60 
and  100  mm.  Two  pressure  transducers  were  installed 
on  the  60-mm-  w  idth  indentor  and  four  pressure  transduc¬ 
ers  on  the  100-mm-width  indentor. 

Support  of  indentor  on  three  load  cells 

The  indentor  plate  was  mounted  on  three  load  cells 
that  were  installed  on  the  supporting  structure  behind 
the  indentor.  The  total  interaction  force  at  the  interface 
was  obtained  by  summing  the  forces  measured  by  each 
load  cell.  The  indentor  widths  used  with  this  setup  were 
50.  100  and  150  mm. 

Mounting  the  indentor  plate  on  three  load  cells 
produced  a  stiff  load-measuring  system  that  gauged  the 
interaction  forces  without  the  influence  of  the  structural 
response.  This  methodology  (Fig.  7b)  to  measure  inter¬ 
action  forces  enabled  us  to  determine  not  only  the 
Figure  5.  Higli-force  module,  looking  from  above.  magnitude  of  the  total  forces  but  also  the  position  of  the 


Figure  6.  Schematic  of  the  experimental  setup. 
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Figure  7.  Schematic  of  indentor  plate. 


resultant  forces.  When  the  resultant  force  is  found  to  be 
in  the  center  of  the  contact  area,  it  is  an  indication  of  a 
symmetrical  distribution  of  ice  pressure  at  the  interface. 
Examples  of  simple  symmetrical  distributions  of  ice 
pressure  would  be  a  uniform  ora  parabolic  distribution 
ora  combination  of  these  two.  However,  if  the  position 
of  the  resultant  force  is  not  at  the  center  of  the  contact 
area,  we  know  that  the  distribution  of  ice  pressure  is 
asymmetrical. 

Instrumentation 

A  44-kN-capacity  (10.000-lb)  load  cell  shown  in 
Figure  6  was  installed  on  the  structural  support  to 
measure  the  total  ice  force.  However,  only  pan  of  the 
total  interaction  force  at  the  interface  was  measured  by 
this  load  cell  (which  was  calibrated  prior  to  the  testing 
program ).  The  capacity  of  each  load  cell  supporting  the 
indentor  plate  was22kN(50()0lb).Thecapacityofeach 
of  the  1 3  pressure  transducers  installed  on  the  indentor 
plate  was  7  MPa  (rated)  and  14  MPa  (maximum). 
Because  the  capacity  of  these  transducers  was  not  high 
enough,  they  failed,  one  after  the  other  during  the  tests 
because  of  overloading. 

The  velocity  of  the  structural  support  was  measured 
using  a  position-velocity  transducer  mounted  on  the 
frame  of  the  high-force  module  and  attached  to  the 
screw-driven  carriage.  Displacement  of  the  indentor 


relative  to  the  ice  sheet  was  measured  using  another 
position-velocity  transducer  that  was  placed  on  the  ice 
sheet  and  connected  to  the  indentor  surface. 

An  Acoustic  Emission  ( AE)  sensor  was  placed  in  the 
ice  sheet  to  monitor  the  acoustic  activity  from  micro¬ 
cracking  during  deformation  of  the  ice.  The  AE 
transducer  was  connected  to  a  signal  conditioner  whose 
output  in  volts  was  proportional  to  the  cumulative  AEs 
(or  counts)  received  by  the  sensor.  It  would  automati¬ 
cally  reset  to  zero  whenever  it  reached  a  maximum 
limit. 

An  event  marker  was  used  to  record  the  initial 
contact  between  the  indentor  and  the  ice  and  the  nuclea- 
tion  of  the  first  radial  macrocrack  in  the  ice  sheet  during 
tests.  The  event  marker  was  connected  to  the  signal 
conditioner,  and  the  researcher  observing  the  tests  cued 
the  signal  conditioner  by  touching  the  event  marker. 

Data  acquisition  system 

A  tlow  chart  of  the  data  acquisition  system  used  to 
monitor  our  tests  is  shown  in  Figure  8.  The  data  acqui¬ 
sition  system  was  controlled  by  a  desk-top  computer 
that  also  provided  input-output  handling  for  data  stor¬ 
age.  All  transducers  were  connected  to  a  signal  condi¬ 
tioner  that  also  provided  the  excitation  voltage  to  each 
transducer.  The  data  were  stored  in  digital  form  on 
floppy  disks. 
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Figure  <S\  Flow  chart  of  rlic  data  acquisition  system. 

Testing  procedures 

When  the  ice  sheet  thickness  reached  that  required 
for  the  planned  test,  the  characteristic  length  of  the  ice 
sheet  was  measured  by  placing  a  deadweight  on  the 
center  of  the  ice  sheet  and  monitoring  its  vertical  dis¬ 
placement  with  respect  to  a  fixed  datum.  A  slot  was  cut 
in  the  ice  sheet  with  a  chain  saw  to  allow  the  indentor. 
which  protruded  below  the  watersurface.  totravel  from 
one  test  location  to  another.  The  ice  was  carefully  cut  to 
avoid  cracking  and  to  make  the  ice  edge  perpendicular 
to  the  indentor  surface.  Indentation  tests  were  per¬ 
formed  on  ice  that  was  visually  free  of  macrocracks. 

The  AE  sensor  was  placed  in  the  ice  sheet  by  drilling 
a  I  -cm-deep  hole  at  a  location  approximately  70  cm  in 
front  of  the  indentor.  A  weight  was  placed  on  it  to  keep 
it  in  firm  contact  with  the  ice  sheet.  A  displacement 
transducer  was  also  placed  on  the  ice  sheet  approxi¬ 
mately  1 .5  m  in  front  of  the  indentor  and  its  string  was 
connected  to  the  indentor  plate.  Indentation  tests  were 
conducted  after  these  preparations. 

Tests  were  conducted  along  several  tracks  in  the  ice 
sheet  as  shown  schematically  in  Figure  9.  First,  a  test 


was  conducted  on  an  undamaged  (with  no  microcracks) 
ice  sheet  for  several  seconds  and  then  stopped  (first 
run).  During  these  tests,  extensive  microcracks  were 
produced  in  front  of  the  indentor.  Then,  anothertest  was 
carried  out  on  this  damaged  ice  in  the  track  created  by 
the  first  run  at  the  same  velocity  (second  run).  The  time 
interval  between  these  two  tests  (or  runs)  was  usually  2 
to  3  minutes.  During  the  first  run.  the  initial  contact  of 
the  indentor  with  the  ice  and  the  nucleation  of  the  first 
radial  macrocrack  were  recorded  manually  with  the 
help  of  the  event  marker. 

After  conducting  two  tests  in  one  track,  we  moved 
the  carriage  to  another  location,  doing  two  more  tests  at 
a  different  velocity,  creating  another  test  track  (see  Fig. 
9).  Up  to  five  test  locations  were  carefully  chosen  to 
avoid  existing  cracks  in  the  ice  sheet,  the  distance 
between  different  test  tracks  being  at  least  3  m.  In  one 
day.  9  to  1 1  tests  were  usually  conducted  on  one  ice 
sheet.  Second  runs  for  tests  at  very  low  velocities  ( I 
mm/s)  were  not  conducted.  Still  photographs  and  video 
recording  of  the  tests  were  taken.  After  the  tests,  crushed 
ice  specimens  were  cut  out  to  observe  and  to  photograph 
the  ice  damaged  by  microcracking. 

Summary  of  test  series 

The  test  period  was  from  19  February  through  18 
April  1 988.  Altogether.  92  tests  were  conducted  in  an  8- 
week  period. 

Tests  were  divided  into  two  test  series.  First,  tests 
were  conducted  in  which  the  pressure  transducers  were 
installed  on  the  indentor  (Fig.  7a).  When  some  of  the 
pressure  transducers  failed  because  of  overloading  during 
the  tests,  the  experimental  setup  was  changed  to  sup¬ 
porting  the  indentor  plate  on  three  load  cells  (Fig.  7b). 
thus  enabling  the  measurement  of  total  force  and  deter¬ 
mining  the  position  of  the  resultant  force.  They  are 
labeled  NSC  test  and  N  N  test,  respectively.  The  test 
conditions  are  listed  in  Table  !  for  NSC  tests  and  in 
Table  2  for  N  N  tests. 


Figure  9.  Test  track  location  in  the  basin. 
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Table  1.  Parameters  and  ice  properties  for  the  NSC 
series  of  tests. 


Dale 
( 19811) 

Test 

number 

huientor 

width 

(mm) 

Ire 

thickness 

(mm) 

Indentor 

velocity 

(mulls) 

Le  (cm)* 

E  (G/ia) 

19  Feb 

89/88 

60 

50 

7.7 

not  measured 

87/86 

60 

42 

6.6 

85/84 

60 

45 

5.4 

83/82 

60 

48 

3.3 

25  Feb 

79 

100 

26 

7.8 

not  measured 

78 

100 

25 

6.6 

77 

100 

25 

4.4 

76 

100 

24 

2.6 

29  Feb 

69/68/67 

100 

57 

7.9 

not  measured 

66 

100 

56 

7.5 

65 

100 

56 

5.5 

64/63/62/61 

100 

56 

5.8 

60 

100 

62 

6.1 

59 

100 

57 

1.8 

7  Mar 

49/48/47 

100 

52 

7.9 

not  measured 

46 

100 

54 

— 

45 

100 

55 

— 

44 

100 

55 

— 

10  Mar 

39/38 

100 

45 

8.0 

not  measured 

37 

100 

43 

8.1 

36 

100 

49 

1.0 

16  Mar 

29 

100 

47 

8.3 

Lc=  116.7 

28 

100 

58 

0.9 

E  =  1.19 

27 

100 

56 

1.5 

26 

100 

56 

1.3 

25 

100 

56 

1.0 

24 

100 

55 

3.1 

23 

100 

53 

4.6 

22 

100 

57 

6.2 

21 

100 

51 

7.9 

18  Mar 

19 

100 

25 

8.0 

£<■  =  65.3 

18 

100 

25 

5.8 

E=  1.25 

17 

too 

25 

3.8 

16 

100 

25 

2.0 

*Lc  =  characteristic  length:  £  =  elastic  modulus. 


TEST  RESULTS  AND  OBSERVATIONS 

In  this  study,  we  observed  microcracking  near  the 
indentor  as  well  as  a  variety  of  failure  modes  during  the 
indentation  tests:  radial  cracking, circumferential  crack¬ 
ing,  cleavage  cracking,  spalling  and  crushing,  and 
buckling.  The  time  when  the  radial  cracks  occurred 
during  a  test  depended  on  the  test  conditions,  i.e.,  ice 
thickness,  indentor  width  and  indentor  velocity. 

Microcracking 

Observations 

Microcracking  was  observed  visually  and  recorded 
by  still  photography  and  video.  Sketches  of  micro- 
cracks  and  some  larger  cracks,  based  on  visual  observa¬ 
tions,  are  shown  in  Figure  10.  Photographs  of  a  test. 


Table  2.  Parameters  and  ice  properties  for  the  N  N 
series  of  tests. 


Dole  Test 
(I0SS)  number 

Indentor 

width 

(mm) 

Ice 

thickness 

(mm) 

Indentor 

velocity 

(minis) 

Lc  (CM)* 

E  (GPa) 

1  Apr 

31/32 

50 

26 

9.3 

£<■  =  65.3 

33/34 

50 

26 

7.0 

£'  =  1.11 

5/36 

50 

27 

4.6 

37/38 

50 

27 

2.2 

39 

50 

27 

1.2 

4  Apr 

41/42 

50 

60 

9.2 

£<•=  146.3 

43/44 

50 

60 

6.8 

£  =  2.27 

45/46 

50 

56 

4.4 

47/48 

50 

57 

2.2 

49 

50 

60 

1.3 

7  Apr 

51 

100 

34 

9.1 

£<■  =  77.3 

52/53 

100 

34 

6.9 

£  =  1 .07 

54 

100 

33 

9.4 

55/56 

100 

33 

4.5 

57/58 

100 

32 

2.3 

59/60 

100 

33 

2.0 

1 1  Apr 

61/62 

100 

53 

9.2 

£<■  =  129.0 

63/64 

100 

56 

6.6 

£  =  1.69 

65/66 

100 

56 

4.6/3.4 

67/68 

100 

54 

2.2 

69/70 

100 

55 

l.l 

701 

100 

57 

9.4 

14  Apr 

71 

150 

30 

4.9 

£<=71.8 

72 

150 

31 

2.9 

£  =  1.05 

73 

150 

30 

0.3 

74 

150 

29 

1.5 

75 

150 

31 

l.l 

76 

150 

33 

1.2 

18  Apr 

81 

150 

54 

3.5 

£<  =  129.0 

E  =  1 .88 


*Lc  =  characteristic  length:  E  =  elastic  modulus. 

which  were  continuously  taken  from  the  initial  contact 
through  the  ice  failure,  are  shown  in  Figure  1 1 . 

At  the  onset  of  loading,  microcracks  formed  in  the 
immediate  vicinity  of  the  indentor  (Fig.  10a).  As  the 
indentor  deformed  the  ice.  more  microcracks  appeared 
in  an  expanding  zone.  The  crack  density  also  appeared 
to  increase  with  the  passage  of  time  after  the  initial 
contact  between  the  indentor  and  the  ice.  As  the  ice 
force  increased  during  microcracking  and  deformation 
of  the  ice.  one  (and  sometimes  two)  big  radial  cracks 
propagated  into  the  ice  sheet  ahead  of  the  indentor.  A 
photograph  of  the  microcracks  and  a  radial  crack  is 
shown  in  Figure  1 2.  With  the  help  of  an  event  marker, 
it  was  possible  to  establish  that  this  radial  crack  formed 
when  the  ice  force  reached  a  magnitude  of  about  70%  of 
the  maximum  force;  there  is  no  indication  in  the  record 
of  any  decrease  in  ice  force  at  the  time  this  crack 
nucleated  and  propagated.  The  length  of  these  radial 
cracks  was  usually  about  3  m  and  they  formed  almost 
perpendicular  to  the  indentor  face  (Fig.  10b). 
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Even  after  the  formation  of  the  radial  crack,  micro- 
cracks  would  keep  developing,  not  only  in  front  of  the 
indentor,  but  also  on  both  sides.  The  ice  in  front  of  the 
indentor  would  turn  milky  with  the  accumulation  of 
microcracks.  We  found  by  observing  the  tests  that  the 
highest  density  of  microcracks  was  near  the  indentor, 
decreasing  with  increasing  distance  from  it  (Fig.  10c). 

When  the  test  velocity  was  high  ( v  >  4  mm/s),  the  ice 
sheet  failed  brittlely,  followed  by  sudden  forward 
movements  of  the  indentor  into  the  ice  sheet  and  into  the 
zone  of  extensive  cracks  (Fig.  1  Od).  When  the  indentor 
velocity  was  low  (v  <  3  mm/s),  the  ice  failed  ductilely, 
as  indicated  by  the  absence  of  sudden  forward  indentor 
motion  and  the  rounded  peaks  in  the  ice  force  record. 

During  tests  at  high  velocity  (v  >  4  mm/s),  we 
observed  spurts  of  microcracking  prior  to  each  subse¬ 
quent  failure.  Microcracks  developed  ahead  of  the 
indentor,  and  then  the  indentor  moved  forward  very 
quickly,  accompanied  by  crushing  or  spalling  of  the 
damaged  (microcracked)  ice  (Fig.  lOe). 


Finnic  /  / .  Progress  of  an  indentation  test 
from  the  initial  contact  through  the  fa  st 
failure. 
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Wonting  it  c  shed 


(i.  d  =100  mm  ami  h  =  54  mm. 


h.  d  =  50  mm  and  h  =  6 0  mm 


Finnic  15.  Acoustic  emission  signals  and  ice  force  records. 


T i me  ( s  ) 


Time  ( s  ) 


a.  Test  N  N  61— d  =  100  mm,  h  =  53 
mm,  v  =  9.2  mmls  . 


b.  Test  N  N  69— d  =  100  mm,  h  =  55 
mm,  v  =  1.1  mmls. 


Figure  16.  Cumulative  AE  count  and  ice  force. 


Microcracks  spread  slower  during  tests  with  a  low 
indentor  velocity  than  they  did  during  tests  at  higher 
velocities.  Figure  13  shows  a  photograph  of  micro- 
cracks  formed  during  a  test,  whereas  Figure  1 4  shows  a 
photograph  of  microcracks  after  a  test.  A  preliminary 
examination  of  thin  sections  from  a  freshly  deformed 
ice  sample  indicated  that  the  microcracks  were  located 
at  the  grain  boundaries.  This  agrees  with  the  observa¬ 
tions  of  Sinha  (1984),  who  proposed  a  model  for  the 
initiation  of  intergranular  cracks  and  their  subsequent 
accumulation. 

Microcracking  and  AE  signals 

Figure  1 5  shows  records  of  the  AE  signals  and  the  ice 
force  versus  time.  As  we  stated  previously  in  the  section 
describing  instrumentation,  the  AE  apparatus  would 
reset  the  cumulative  count  (or  the  signal)  automatically 
to  zero  when  the  count  reached  the  saturation  level  of 
the  apparatus.  The  i ntense  AE  act  i  v ity  i s  reflect ed  i n  the 
records  by  the  high  frequency  of  resetting.  The  increase 
in  ice  force  correlates  very  well  with  the  high  AE  count, 
which  can  be  directly  attributed  to  microcracking  as 
observed  visually.  For  the  ductile  failure  at  low  velocity 
(v  <  3  mm/s),  an  AE  signal  was  recorded  after  the  time 
when  the  maximum  ice  force  occurred. 

Figure  1 6  shows  the  cumulative  AE  count  and  the  ice 
force  with  respect  to  the  time  for  the  indentor  velocities 
of9.2and  1.1  mm/s.Theuppercurvesinthe  AEfigures 
are  the  cumulative  count,  which  has  been  nondimen- 
sionalized  with  respect  to  its  maximum  value.  After  the 
buildup  in  the  ice  force,  the  cumulative  AEcount  started 
to  increase,  and  the  ice  force  and  the  cumulative  AE 
count  increased  together.  These  findings  show  that  the 
microcracks  are  generated  and  accumulated  in  the  ice 
by  the  indentor  loading.  In  the  test  having  high  indentor 


velocity  (9.2  mm/s),  the  cumulative  AE  count  increased 
until  brittle  failure.  At  the  low  indentor  velocity  ( 1 . 1 
mm/s),  however,  the  cumulative  AE  count  increased 
with  increasing  ice  force.  Even  after  the  time  when  the 
peak  force  was  recorded,  this  count  continued  to  in¬ 
crease  because  of  the  ductile  behavior  of  the  ice. 

Failure  modes 

The  failure  modes  observed  in  the  test  series  are 
illustrated  schematically  in  Figure  17. 

Radial  crack 

As  mentioned  in  the  previous  section,  one  (or  some¬ 
times  two)  radial  cracks  were  generated  after  a  certain 
time  (between  1  and  5  seconds,  which  depended  on 
indentor  velocity)  following  the  initial  contact  between 
the  indentor  and  the  ice.  The  length  of  these  radial 
cracks  was  usually  about  3  m  and  they  formed  almost 
perpendicular  to  the  indentor  face  (Fig.  17a). 

Circumferential  crack 

Many  circumferential  cracks  were  generated  in  front 
of  the  indentor  by  the  first  brittle  failure  (Fig.  17b). 
Usually,  five  to  ten  circumferential  cracks  developed  in 
front  of  the  indentor.  No  circumferential  cracks  were 
generated  when  the  ice  failed  ductilely. 

Cleavage  crack 

Cleavage  cracks  (in-plane  horizontal  cracks)  were 
sometimes  observed  in  front  of  the  indentor.  These 
cracks  appeared  afterthe  microcracks  formed  and  grew 
gradually  in  size  as  the  indentor  loading  increased  (Fig. 
1 7c).  In  Figure  1 8.  the  dense,  milky  semicircular  shape 
in  front  of  the  indentor  is  a  cleavage  crack.  Figure  19 
shows  a  photograph  of  both  a  plan  view  and  a  cross 
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Figure  17.  Global  failure  modes  (after  Sanderson  1988). 


section  of  the  ice  cut  after  the  test  shown  in  Figure  1 8.  lead  to  tensile  forces  at  the  crack  tip.  More  research  is 

A  theoretical  model  for  the  propagation  of  cleavage  needed  to  understand  cleavage  crack  propagation, 

cracks  was  presented  by  Kendall  (1978).  A  geometry 
useful  in  the  analysis  of  cleavage  propagation  is  shown  Spalling 

in  Figure  20.  The  compressive  force  P  is  conceptually  During  tests  at  intermediate  to  high  indentation  ve- 

di  vided  into  two  forces  P/2  acting  at  points  a  distance  III  locities  (r  >  6  mm/s).  we  observed  local  spalling  failure 

2  from  the  top  and  bottom  surfaces.  These  two  forces  during  ice  crushing.  The  sudden  forward  motion  of  the 


Figure  19.  Ice  horizontally  separated  by  the  in-plane  cleavage  crack: 
the  plan  view  and  the  cross  section  of  the  ice  (arrow  indicates 
indentor  path). 


a.  Plan  view.  b.  Cross-sectional  view for  the  analysis 

of  propagation  of  the  cleavage  crack. 

Figure  20.  Geometry  of  an  in-plane  cleavage  crack  (after  Wierzbicki  1985). 


Figure  21 .  Buckled  ice  sheet  with  a  number  of  circumferential  cracks. 


indentor  caused  fragments  of  ice  to  break  away  ahead  of 
it.  Ice,  damaged  by  niicrocracks.  expanded  near  the 
indentor.  As  the  ice  moved  in  the  vertical  direction,  the 
indentor  moved  forward  suddenly,  with  the  ice  frag¬ 
ments  spalling  up  and  down  (Fig.  17d). 

Buckling 

In  some  tests,  ice  buckled,  especially  when  indentors 
were  pushed  against  thin  ice  at  low  velocities  ( v  <  3  mm/ 
s).  Tests  N  N  7 1  to  N  N  76  (c/  =  150  mm)  all  resulted  in 
buckling  failure.  At  the  onset  of  the  loading,  micro¬ 
cracks  were  generated  in  the  same  manner  as  stated 
earlier.  But  soon  after  the  nucleation  of  the  microcracks, 
the  ice  sheet  started  to  bend  downward  and  eventually 
buckled  forming  several  circumferential  cracks  (Fig. 
21). 


DISCUSSION 

To  examine  the  factors  influencing  the  ice  forces  on 
flat,  vertical  indentors,  we  looked  at  the  following 
topics:  velocity  effects,  energy  analysis,  ice  pressure 
and  contact  area,  position  of  resultant  force,  first  peak 
force  versus  subsequent  peak  force,  frequency  of  ice 
force  failure,  and  buckling  failure. 

The  velocities  indicated  in  figures  and  tables  are 
those  of  the  average  rate  of  indentor  displacement 


relative  to  the  ice  sheet.  Test  results  and  analyses  are 
summarized  in  Table  3. 

Velocity  effects 

Because  ice  is  a  viscoelastic  material,  its  properties 
depend  on  the  magnitude  of  the  stress  it  is  underand  the 
rate  of  stress  application.  The  ice  may  fail  by  brittle 
fracture  or  by  ductile  creep,  depending  on  the  applied 
stress  rate.  The  range  of  relative  velocity  in  the  present 
test  series  was  between  1  and  9  mrn/s. 

At  the  instant  of  ice  failure,  especially  during  tests  at 
intermediate  to  high  indentor  velocity  (v  >  4  mrn/s ).  a 
sudden  unloading  occurred  with  a  loud  noise  because 
the  indentor  would  spring  back  to  its  original  position. 
This  can  be  seen  in  the  displacement  records  of  the 
carriage  and  the  indentor  in  Figure  22,  where  force  and 
displacement  records  are  shown  for  tests  at  velocities 
ranging  from  3.1  to  7.9  mm/s.  Figure  23  defines  the 
deformation  of  the  structure:  the  difference  between  the 
displacement  of  the  carriage  and  that  of  the  indentor. 
We  observed,  from  Figure  22.  that  the  indentor  dis¬ 
placement  record  was  not  always  linear  with  time.  The 
nonlinear  displacement  record  of  the  indentor  at  a 
velocity  of  3.1  mm/s  in  Figure  22  shows  that  the 
indentor  catches  up  with  the  carriage  faster  as  the  ice 
resistance  decreases  as  a  result  of  creep  failure. 

The  effects  of  velocity  on  the  nature  of  failure  can  be 
seen  in  the  records  of  ice  force  versus  time  and  in  the 
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records  of  ice  force  versus  indentor  displacement  shown 
in  Figures  24  and  25  respectively.  The  sudden  unload¬ 
ing  of  the  force  on  the  indentor  is  an  indication  of  brittle 
failure,  whereas  the  ductile  failure  results  in  gradually 
decreasing  force,  but  not  all  the  way  to  zero.  Typical 
force  versus  time  plots  for  the  brittle  and  the  ductile 
failures  of  ice  are  shown  in  Figure  26. 

At  high  velocities  (v  >  4  mm/s),  after  the  initial 
contact  of  the  indentor  with  the  ice,  the  ice  force 
increases  monotonically  to  a  maximum  value,  as  shown 


in  Figures  24  and  25.  Microcracks  nucleate  and  accu¬ 
mulate  in  the  ice  sheet  during  this  period.  Brittle  ice 
behavior  results  in  sudden  unloading  of  the  indentor  at 
the  instant  of  the  ice  failure  and  forward  indentor 
movement,  with  crushing  or  spalling  of  the  micro- 
cracked  ice  in  front  of  the  indentor.  Because  of  brittle 
failure,  a  sawtooth  type  of  force-time  plot  is  recorded 
for  high  velocity  tests. 

At  low  velocities  (v<  3  mm/s),  the  data  in  Figures  24 
and  25  show  that  the  ice  forces  gradually  increase  to  a 


Table  3.  Summary  of  test  parameters  and  results. 


Ice 

Elastic 

Indentor 

R  = 

thickness,  h 

modulus.  E 

Width.  4  Velocity,  v 

F^* 

p  ** 

F.J 

(mm) 

(GPa) 

(mm)  (mm/s) 

(kN) 

(MPa) 

(kN) 

fh/v 

NSC 

21 

51 

1.19 

100 

7.9 

70.0 

13.7 

31.8 

0.45 

22 

57 

1.19 

100 

6.2 

67.1 

1 1.8 

29.3 

0.44 

0.54 

4.96 

23 

53 

1.19 

100 

4.6 

63.0 

1 1.9 

24 

55 

1.19 

100 

3.1 

63.0 

11.5 

25 

56 

1.19 

100 

1.0 

56.0 

10.0 

NN 

31/32 

26 

1.11 

50 

9.3 

17.0 

13.1 

6.4 

0.38 

33/34 

26 

1.1 1 

50 

7.0 

15.6 

12.0 

5.8 

0.37 

35/36 

27 

1.1 1 

50 

4.6 

15.3 

11.3 

37/38 

27 

1.11 

50 

2.2 

13.7 

10.1 

39 

27 

l.ll 

50 

1.2 

12.7 

9.4 

41/42 

60 

2.27 

50 

9.2 

39.7 

13.2 

20.9 

0.53 

0.75 

4.89 

43/44 

60 

2.27 

50 

6.8 

38.4 

12.8 

20.5 

0.53 

0.59 

5.21 

45/46 

56 

2.27 

50 

4.4 

35.7 

12.8 

47/48 

57 

2.27 

50 

2.2 

33.0 

1 1.6 

49 

60 

2.27 

50 

1.3 

31.5 

10.5 

51 

34 

1.07 

100 

9.1 

31.9 

9.4 

13.0 

0.41 

1.52 

5.68 

52/53 

34 

1.07 

100 

6.9 

28.4 

8.4 

1 1.0 

0.39 

0.86 

4.24 

54 

33 

1.07 

100 

9.4 

28.7 

8.7 

13.0 

0.45 

1.38 

4.84 

55/56 

33 

1.07 

100 

4.5 

28.7 

8.7 

57/58 

32 

1.07 

100 

2.3 

22.0 

6.9 

59/60 

33 

1.07 

100 

2.0 

24.2 

7.3 

61/62 

53 

1.69 

100 

9.2 

57.2 

10.8 

18.0 

0.31 

0.94 

5.42 

63/64 

56 

1.69 

100 

6.6 

54.6 

10.3 

23.8 

0.44 

0.57 

4.84 

65/66 

56 

1.69 

100 

4.0 

52.6 

9.4 

23.3 

0.44 

0.38 

5.32 

67/68 

54 

1.69 

100 

2.2 

46.9 

8.7 

69/70 

55 

1.69 

100 

1.1 

42.5 

7.7 

701 

57 

1.69 

100 

9.4 

52.1 

9.1 

28.7 

0.55 

0.72 

4.37 

71 

30 

1.05 

150 

4.9 

39.3 

8.7 

72 

31 

1.05 

150 

2.9 

31.0 

6.7 

73 

30 

1.05 

150 

0.3 

25.0 

5.6 

74 

29 

1.05 

150 

1.5 

33.2 

7.6 

75 

31 

1.05 

150 

l.l 

28.3 

6.1 

76 

33 

1.05 

150 

1.2 

30.5 

6.2 

*  Maximum  force. 

**  Effective  pressure, 
t  Maximum  force  of  the  first  peak, 
tt  Frequency  of  ice  failure. 
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F  inure  22.  Plots  of  ice  force  and  displacement  of  the  carriage 
and  the  indentor  with  respect  to  rime  for  the  tests  with  d  =  100 
mm  and  h  =  54  mm  at  different  velocities. 


maximum  value  after  the  initial  contact  and  that  the 
ductile  deformation  of  ice  results  in  gradual  unloading 
of  the  indentor  and  smooth  variations  in  the  ice  force 
records.  Ductile  failure  can  be  recognized  by  rounded 
peaks  attributable  to  softening  of  the  material  at  high 
stress  levels. 

Sharp  peaks  in  force  versus  time  plots  become 
rounded  with  decreasing  velocity.  It  should,  however, 
be  noted  that  even  at  high  velocity  (v  >  4  mm/s),  the 
monotonic  increase  in  the  ice  force  becomes  slightly 
rounded  before  the  ice  failure  takes  place. 

Further,  the  plots  of  ice  force  versus  the  indentor 
displacement  at  different  indentor  velocities  (Fig.  27) 
present  ice  force  records  up  to  the  first  few  peaks  to 
show  the  effect  of  velocity. 

From  these  records,  it  appears  that  the  indentor 
displacements  corresponding  to  the  peak  ice  force  are 
almost  the  same  even  though  the  indentor  velocities  are 
different  (the  peaks  in  the  records  for  ductile  failure  are 
a  little  bit  to  the  right  of  those  for  brittle  failure).  Also, 
even  after  the  ice  failed  brittlely  (i.e„  tests  N  N  35, 43, 
45, 52  and  6 1 ),  the  indentors  are  still  loaded  for  a  certain 
distance  after  the  peak  force  is  recorded. 

After  the  sudden  unloading  of  the  indentor,  it  moves 
forward  until  it  again  contacts  the  ice.  Figure  27  shows 
that  the  second  loadings  start  almost  at  the  same  posi¬ 
tion  in  the  ice  sheet  for  different  indentor  velocities.  The 


size  of  the  damage  zone  caused  by  the  first  peak  loading 
is  about  the  same  even  when  the  indentor  velocities  are 
different. 

The  maximum  effective  pressure,  defined  as  the 
maximum  ice  force  divided  by  the  contact  area,  is 
plotted  with  respect  to  the  indentation  rate  v/d  (i  = 
indentor  velocity,  d  =  indentor  width)  in  Figure  28.  We 
can  see  a  trend  of  the  maximum  effective  pressure 
slightly  increasing  with  the  increase  in  indentation  rate 
vld,  which  is  in  the  range  of  0.0 1  and  0. 1 8  1/s. 


Figure  23.  Sketch  explaining  the  dis¬ 
placement  records  shown  in  Figure 
22:  a — displacement  of  the  carriage: 
b — displacement  of  the  indentor  with 
respect  to  the  ice  sheet:  c — deforma¬ 
tion  of  the  structure. 


indcntor  displacement  records  at  different  velocities. 


Force  CkN)  Force  (kN) 


a.  Brittle  failure;  saw  tooth 
type  records  with  force  going 
to  zero  after  each  failure. 


h.  Ductile  failure;  rounded  ice 
force  curve  with  force  nor  going 
to  zero. 


Figure  26.  Typical  ice  force  records. 
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Figure  27.  Comparison  of  ice force  versus  indentor  displacement  records  up  to  the  first  few  peaks  at  different  velocities. 
The  upper  five  plots  show  separately  the  records  for  different  velocities;  the  bottom  plot  shows  the  records  together. 
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In  Figure  29,  plots  of  maximum  effective  pressure 
versus  indentation  rate  obtained  from  the  present  study 
are  compared  with  those  from  other  studies  on  labora¬ 
tory  ice  (Hirayama  et  al.  1974,  Frederking  and  Gold 
1975,  Michel  and  Toussaint  1977,  Timco  1987).  We 
obtained  an  average  value  for  maximum  effective  pres¬ 
sure  of  about  10  MPa,  which  is  consistent  with  the 
results  of  the  previous  studies. 

Energy  analysis 

Energy  supplied  to  the  carriage  by  the  drive  screw  is 
partly  stored  in  the  structure  and  partly  dissipated  in  the 
ice.  Energy  stored  in  the  structure  results  in  the  defor¬ 
mation  of  the  structure,  and  the  energy  dissipated  in  the 
ice  causes  deformation  and  microcracking  in  the  ice. 
The  force  generated  at  the  interface  is  the  same  as  that 
which  deforms  the  ice,  deflects  the  structure  and  moves 
the  carriage.  Because  the  displacements  of  the  carriage 
and  the  indentor  were  measured  separately  with  respect 
to  a  fixed  datum,  the  energies  supplied  to  the  carriage 
and  dissipated  in  the  ice  can  be  computed  independ¬ 
ently.  The  energy  stored  in  the  structure  is  computed  as 
the  difference  between  the  energy  supplied  to  the  car¬ 
riage  and  the  energy  dissipated  in  the  ice. 

Figure  30  shows  measured  ice  force  and  displace¬ 


ments  data  and  computed  energy  transfers  for  indentor 
velocities  of  9.4  and  1 .3  mm/s  respectively.  It  shows 
plots  of  the  following:  (from  the  top)  the  ice  force;  the 
displacements  of  the  carriage  and  the  indentor,  the 
energies  supplied  to  the  carriage,  dissipated  in  the  ice 
and  stored  in  the  structure;  and  the  rates  of  energy 
supplied  to  the  carriage,  stored  in  the  structure  and 
dissipated  in  the  ice  with  respect  to  the  time.  Energy 
supplied  to  the  carriage  was  calculated  by  integrating 
the  product  of  the  ice  force  and  the  incremental  carriage 
displacement.  Energy  dissipated  in  the  ice  sheet  was 
calculated  by  integrating  the  product  of  the  ice  force  and 
the  incremental  indentor  displacement.  Energy  stored 
in  the  structure  was  calculated  by  subtracting  the  energy 
dissipated  in  the  ice  from  the  energy  supplied  to  the 
carriage.  We  calculated  energy  rates  by  multiplying  the 
force  with  the  corresponding  velocities,  which  in  turn 
were  obtained  by  fitting  a  polynomial  through  the 
displacement  data  by  the  method  of  least-squares  and 
then  differentiating  the  fitted  polynomial  with  respect  to 
time. 

From  the  initial  contact  to  a  peak  force,  a  portion  of 
the  energy  from  the  carriage  is  stored  in  the  structure  as 
the  structure  deforms  as  a  result  of  the  interaction  force. 
During  a  test  at  high  indentor  velocity  (e.g..  9.4  mm/s. 
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Energy  Rite  (H)  Energy  (J)  Displ.(mm)  Force  (kN) 


a.  Test  N  N  41 ,  v  =  9.4  minis.  b.  Test  N  N  49,  v  =  1.3  minis. 

Figure  30.  Energy  and  energy-rate  plots  (d  =  50  mm  and  h  =  60  mm). 


Fig.  30a),  the  energies  stored  in  the  structure  and  dissi¬ 
pated  in  the  ice  increase  with  the  elapsed  time  until  the 
force  reaches  a  maximum  value.  At  this  time,  the  energy 
stored  in  the  structure  is  greater  than  that  dissipated  in 
the  ice.  This  means  that  the  deformation  of  the  structure 
is  larger  than  the  displacement  of  the  indentor  relative 
to  the  ice  sheet  until  the  ice  fails. 

During  a  test  at  low  velocity  (e.g.,  1.3  mm/s.  Fig. 
30b),  the  energy  in  the  structure  and  in  the  ice  increases 
with  the  elapsed  time,  but  the  energy  dissipated  in  the 
ice  is  greater  than  that  stored  in  the  structure.  This  means 
that  the  displacement  of  the  indentor  relative  to  the  ice 
sheet  is  greater  than  the  deformation  of  the  structure 
during  the  First  loading. 

The  rates  of  energy  supplied  to  the  carriage,  stored  in 
the  indentor  structure  and  dissipated  in  the  ice  are 
shown  at  the  bottom  of  Figure  30.  For  a  test  with  high 
indentor  velocity  (v  =  9.4  mm's),  the  rate  of  energy 
dissipated  in  the  ice  is  low  and  gradually  increases  with 
the  elapsed  time;  the  curve  becomes  steep  at  the  time  of 
the  ice  failure.  Fora  test  with  low  indentor  velocity  (v  = 

1 .3  mm/s),  the  rate  of  energy  dissipated  in  the  ice  is 
much  greater  than  that  stored  in  the  structure,  and  it  is 
even  greater  than  that  supplied  to  the  carriage  after  the 


force  peaks  and  the  structure  unloads.  The  negative 
values  for  the  rate  of  energy  stored  in  the  structure 
means  that  the  structure  gradually  swings  back  to  its 
original  position,  thereby  losing  its  stored  energy. 

Appendix  A  presents  20  plots  showing  measured 
data  and  computed  energy  transfers  (Fig.  A92-A  111). 

Figure  31  shows  plots  of  time  versus  dissipation 
energy  in  the  ice  up  to  the  first  peak  in  ice  force  for 
different  indentorvelocities.  In  each  test  series,  itshould 
be  noted,  once  again,  that  the  energy  used  to  produce  the 
peak  force  is  almost  the  same,  even  when  the  indentor 
velocity  is  different.  The  reason  for  this  result  can  be 
seen  in  the  force-displacement  records  shown  in  Figure 
27 .  With  the  lower  indentor  velocity,  the  the  peak  forces 
are  slightly  smaller,  and  these  occur  slightly  to  the  left, 
resulting  in  a  constant  area  under  the  force-displacement 
plots.  In  other  words,  although  the  peak  force  and  the 
displacement  of  the  indentor  are  different  at  different 
indentor  velocities,  the  energy  dissipated  in  the  ice  up  to 
the  first  peak  in  force  is  almost  the  same  if  indentor 
width,  ice  thickness  and  structural  stiffness  are  the 
same. 

Rates  of  energy  dissipated  in  the  ice  up  to  the  first 
peak  in  force  fordifferent  indentor  velocities  are  plotted 
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Energy  Rate  (W)  Energy  Rate  <H)  Energy  (J) 


a.  d  =  50  mm.  b.  d  =  100  mm. 

Figure  3 1 .  Energy  dissipated  in  the  ice  versus  time  far  different  indentor  velocities  (number  at  the  top  of  each  plot). 


in  Figure  32  for  each  test  series.  Each  plot  shows  that  the 
rate  of  energy  dissipated  in  the  ice  increases  with  the 
elapsed  time.  For  tests  at  high  indentor  velocity  (v  >  4 
mm/s),  the  rate  increases  rapidly  with  the  elapsed  time. 
The  rates  are  much  lower  for  tests  at  low  indentor 
velocity  (v<  3  mm/s)  and  attain  a  steady-state  value  as 
opposed  to  the  steep  rise  observed  for  high-velocity 
tests.  This  high  rate  of  energy  dissipation  at  high  inden¬ 
tor  velocity  leads  to  instability,  which  is  commonly 
called  brittle  failure  of  ice.  The  low  rate  of  energy 
dissipation  in  the  ice  at  low  indentor  velocities  leads  to 
gradual  unloading  of  the  indentor,  which  is  commonly 
called  the  ductile  failure  of  the  ice. 


Icc  Th t ckn«**-2?mm  Ic«  Th  ( ckn«s*-60mm 


Zca  Th I ckn***-33mm  Ice  Th I cknet «“55mm 


b.  d  =  100  mm. 


Figure  32.  Energy  rate  dissipated  in  the  ice  versus  time 
for  different  indentor  velocities  (number  at  the  top  of 
each  plot). 


Ice  pressure  and  contact  area 

From  a  wide  range  of  earlier  experiments  on  the 
indentation  of  ice,  it  is  known  that  the  effective  pressure 
decreases  with  increasing  aspect  ratio  or  contact  area. 
Ice  in  the  laboratory  fails  at  10-20  MPa  and  yet  at  large 
scale  in  the  Field  it  appears  to  fail  at  about  1  MPa.  It  is 
clear  that,  for  some  reason,  failure  pressure  depends  on 
contact  area  (Sanderson  1986). 

First,  the  effective  pressures  obtained  from  this  study 
were  examined  to  see  whether  they  followed  any  trend. 
Second,  the  pressure  transducer  records  were  analyzed 
to  see  if  there  was  non-simultaneous  failure  across  the 
contact  area. 

Contact  area  effects 

Figure  33  shows  the  plots  of  the  maximum  effective 
pressure  p  (=  F/dlu  where  F  =  total  force,  d  =  indentor 
width  and//  =  ice  thickness)  as  a  function  of  contact  area 
(d  x  li)  from  the  present  study  for  indentor  widths  of  d 
=  50,  100  and  1 50  mm.  The  test  data  do  not  show  any 
trend  of  decreasing  ice  pressure  with  increasing  contact 
area,  probably  because  of  the  small  variations  in  contact 


Figure  33.  Maximum  peak  pressure  versus  contact  area 
in  the  present  study.  The  numbers  beside  each  point  are 
the  indentor  velocity  in  millimeters  per  second. 
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Figure  35.  Maximum  peak  pressure  versus  aspect  ratio  for  all 
published  tests  that  used  a  flat  indentor  at  a  high  loading  rate 
(after  Timco  1987). 
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Figure  36.  Measured  pressure  records  along  with  force 
records,  d  =  100  mot  and  h  -54  mm  at  different  velocities. 


Figure  36  shows  the  pressure  transducer  records  along 
with  the  ice  force  records  for  tests  in  which  the  indentor 
width  was  100  mm  and  ice  thickness  54  mm.  The 
velocity  range  in  this  figure  is  between  3. 1  and  7.9  mm/ 
s.  Nearly  uniform  pressure  distribution  can  be  seen 
during  initial  contact  of  the  indentor  with  the  edge  of  an 
ice  sheet.  In  a  few  tests,  the  interfacial  pressure  in¬ 
creased  to  a  certain  value  (approximately  1 0  MPa  in  Fig. 
36)  soon  after  the  indentor  contacted  the  ice  and  re¬ 
mained  more  or  less  constant  while  the  ice  force  in¬ 
creased  gradually  until  the  ice  failed. 

A  number  of  pressure  transducer  records  are  in¬ 
cluded  in  NSC  test  sheets  in  Appendix  A.  Although 
some  of  the  records  show  non-simultaneous  crushing 
behavior,  we  believe  that  pressure  data  from  more 
pressure  transducers  will  be  needed  to  ascertain  the  non- 
simultaneous  failure  of  ice  from  the  pressure  distribu¬ 
tion. 

Position  of  resultant  force 

As  mentioned  earlier,  the  indentor  plate  was  mounted 
on  three  load  cells  after  the  change  in  the  experimental 
setup.  Mounting  the  indentor  plate  directly  on  the  load 
cells  is  a  better  way  to  measure  the  interaction  force  than 
that  by  instrumenting  the  structural  support.  This  setup 
allowed  us  to  obtain  not  only  the  total  ice  force  by 
summing  the  three  ice  forces  obtained  from  the  load 
cells  but  also  the  position  of  the  resultant  force  in  the 
contact  area. 


The  positions  of  the  resultant  force  were  calculated 
by  the  following  equations 


r*  ■ 

2  (F1+F2+F3) 


y  = - El - H 

f]+f2+f3 


where  F{,  F2 and  F}  =  ice  force  measured  by  load 
cells  Ly,  L2  and  L} 

B  =  horizontal  distance  be¬ 
tween  load  cells  and  L, 

H  =  vertical  distance  between  the 
Z^-L,  line  andZ.v 

The  calculated  positions  of  the  resultant  force  within 
the  contact  area  are  plotted  in  Figure  37  for  the  total 
force  exceeding  certain  threshold  values:  5  kN  in  Fig. 
37a,  15  kN  in  Fig.  37b.  10  kN  in  Fig.  37c  and  30  kN  in 
Fig.  37d. 
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b.  d  =  50  mm,  h  =  60  mm for  resultant  force 
>  15  kN. 
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Figure  37.  Resultant  force  position  as  determined  from  force  measurements  by  three  load  cells  (resultant  forces  are 
shown  for  each  test  on  the  upper  graph). 
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Figure  37  (cont'cl).  Resultant  force  position  as  determined  from  force  measurements  by  three  load  cells  (resultant 
forces  are  shown  for  each  test  on  the  upper  graph). 
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Figure  38.  ice  force  records  at  different  velocities  of  fust  run  and  second  run  in  the  same  track. 


By  monitoring  the  point  of  action  of  the  resultant 
force,  a  few  inferences  can  be  made  about  the  pressure 
distribution.  A  symmetrical  pressure  distribution  would 
be  indicated  by  the  point  of  action  of  the  resultant  force 
being  in  the  center  of  the  contact  area;  otherwise,  we 


would  expect  it  to  be  in  some  other  area  of  the  contact 
zone.  In  the  situation  of  a  non-simultaneous  ice  failure, 
the  position  of  the  resultant  force  is  expected  to  move 
from  one  part  of  the  contact  zone  to  another  during 
continuous  crushing  of  ice. 
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Figure  39.  Ratio  of  the  mean  subsequent 
peak  force  to  first  peak  force  plotted  versus 
d/h,  dxh  and  v/h. 


The  plots  of  resultant  force  positions  in  Figure  37 
show  that  the  initial  contact  point  between  the  indentor 
and  the  ice  edge  may  not  have  been  in  the  center  of  the 
contact  area,  but,  as  the  load  increased,  the  position  of 
the  resultant  force  moved  to  the  center  of  the  contact 
area.  From  our  observations  of  the  computer  graphics 
during  data  analysis,  weconsistently  saw  the  position  of 
the  resultant  force  move  to  the  center  of  the  contact  zone 
as  the  ice  force  increased.  These  records  indicate  that, 
when  the  interaction  forces  are  greater  than  the  thresh¬ 
old  force  level,  a  symmetrical  pressure  distribution 


exists  at  the  ice-indentor  interface.  It  appears  that  the 
position  of  the  resultant  force  at  high  velocity  (v  >  4  mm/ 
s)  is  slightly  scattered  because  of  the  cyclic  brittle 
failure.  At  low  velocities  (v  <  3  mm/s),  the  position  is 
laterally  centered  but  varies  somewhat  vertically.  This 
may  be  attributed  to  the  gradual  loading  and  the  ductile 
behavior  of  the  ice,  generating  a  symmetrical  pressure 
distribution. 

It  should  be  noted  that  the  uniformity  of  the  ice 
pressure  may  be  attributed  to  the  small  contact  area  in 
this  testing  program.  Further  studies  employing  larger 
contact  areas  should  be  conducted  to  investigate  non- 
simultaneous  failure  of  ice  within  the  contact  area. 

First  peak  force  versus  subsequent  peak  force 

We  stated  previously  that  eight  to  ten  tests  were 
conducted  in  pairs  in  four  to  five  tracks  in  an  ice  sheet. 
Two  tests  were  conducted  in  one  test  track  at  the  same 
velocity  to  observe  the  magnitude  of  the  ice  force 
resulting  from  interaction  with  “undamaged”  and 
“damaged”  ice:  the  “first  run”  into  undamaged  ice  and 
the  "second  run”  into  damaged  ice. 

The  ice  force  records  are  shown  in  Figure  38.  Tests 
N  N  31  and  N  N  32  (Fig.  38a)  were  conducted  in  the 
same  test  track,  as  were  N  N  33  and  N  N  34  and  so  on. 
The  first  peak  force  in  the  first  run  was  always  the 
maximum  force  throughout  each  test.  The  first  peak  was 
generated  by  the  interaction  between  the  indentor  and 
the  undamaged  (with  no  microcracks)  ice  sheet,  and 
subsequent  peak  forces  were  generated  by  indentation 
into  the  damaged  ice.  In  the  second  run,  the  peak  forces 
were  not  as  large  as  the  first  peak  force  in  the  first  run. 
These  lower  peak  forces  in  the  second  run  can  be 
attributed  to  the  previous  damage  and  microcracking 
present  in  ice  caused  during  the  first  run. 

The  ratio  of  the  average  subsequent  peak  force  to  the 
first  peak  force  (F  /F^J  is  plotted  with  respect  to  the 
aspect  ratio  d/h  in  Figure  39a,  the  contact  area  d  x  li  in 
Figure  39b  and  the  velocity-to-thickness  ratio  v/h  in 
Figure  39c.  In  each  of  three  plots,  the  ratio  F  JF  is 
in  the  range  of  0.3  to  0.6,  and  these  plots  show  no  in¬ 
creasing  or  decreasing  trend. 

Frequency  of  ice  force  failure 

Analysis  of  the  structural  vibration  caused  by  re¬ 
peated  ice  failure  is  an  important  part  of  the  design 
process.  Peyton  (1968)  and  Blenkam  (1970)  reported 
on  the  vibrations  of  structures  caused  by  ice  movement 
in  the  Cook  Inlet,  Alaska.  These  structures  were  multi- 
legged  drilling  platforms  that  were  instrumented  to 
measure  ice  forces  and  structure  response.  Jefferies  and 
Wright  (1988)  have  reported  on  the  ice-induced  vibra¬ 
tions  of  Molikpaq,  a  large  structure  (116x116  m)  that 
was  placed  on  a  14-m-deep  berm  at  the  Amauligak 
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Figure  40.  Ice  failure  frequency  plotted 
versus  v,  v/h  and  d/h. 


location.  Maattanen  ( 1 975)  conducted  extensive  meas¬ 
urement  of  the  ice-induced  vibrations  of  lighthouses  in 
the  Gulf  of  Bothnia. 

As  can  be  seen  from  the  data  in  Figures  24  and  25,  the 
interaction  force  between  the  indentor  and  the  ice  in¬ 
creases  with  time  (or  displacement)  until  the  ice  fails 
and  is  not  able  to  resist  the  indentor.  As  has  been 
demonstrated,  depending  on  the  velocity  of  the  inden¬ 
tor,  the  failure  in  the  ice  sheet  is  either  brittle  or  ductile. 
For  brittle  failure  of  ice  at  high  indentor  velocity  (v  >  4 
mm/s),  there  is  sudden  unloading  of  the  indentor;  thus, 
a  brittle  failure  event  in  ice  is  associated  with  the 
increase  of  ice  force  and  a  sudden  unloading. 


After  the  ice  failure,  the  indentor  moves  forward, 
extruding  the  crushed  ice  in  front  of  it.  The  ice  force 
does  not  rise  until  the  indentor  again  contacts  intact  ice. 
The  forward  distance  moved  by  the  indentor  during 
each  failure  event  can  be  characterized  as  a  damage 
zone  caused  by  crushing.  The  data  in  Figure  25  indicate 
that  the  greater  the  ice  force  recorded,  the  farther  the 
indentor  moved  forward.  In  other  words,  the  size  of 
damage  zone  during  successive  failure  events  depends 
on  the  magnitude  of  the  ice  force.  But  an  average 
damage  zone  can  be  obtained  by  counting  the  number  of 
peaks  in  the  force-displacement  record  and  dividing  it 
by  the  total  indentor  displacement.  An  alternate  proce¬ 
dure  to  arrive  at  the  average  damage  zone  is  given 
below. 

Frequencies  of  ice  failure  (/)  in  the  brittle  range, 
listed  in  Table  3,  were  calculated  by  counting  the 
number  of  peaks  divided  by  the  elapsed  time.  These 
frequencies  are  plotted  with  respect  to  the  indentor 
velocity  in  Figure  40a  for  all  tests  with  brittle,  repeated 
failures,  showing  a  trend  of  increasing  frequency  with 
increasing  velocity.  However,  the  scatter  in  the  plots 
increases  with  velocity.  When  the  frequency  is  plotted 
with  respect  to  the  velocity-to-thickness  ratio  v/h  (Fig. 
40b),  the  variation  in  the  data  points  reduces  and  the 
frequency  increases  linearly  with  increasing  v/h.  These 
trends  are  similar  to  those  obtained  by  Sodhi  and  Morris 
(1984). 

A  plot  of  frequency  versus  the  aspect  ratio  d/h  is 
shown  in  Figure  40c.  Although  it  seems  that  the  fre¬ 
quency  increases  with  increasing  the  aspect  ratio,  it  is 
hard  to  ascertain  this  trend  because  of  the  small  vari¬ 
ations  in  the  aspect  ratio. 

The  nondimensional  variable  v/fh  is  the  ratio  of  the 
damage  zone  to  ice  thickness,  where  v  is  the  indentor 
velocity,  h  the  ice  thickness  and /the  frequency  of  ice 
failure.  This  number  quantifies  the  distance  that  the 
indentor  moves  forward  in  each  loading  cycle  in  terms 
of  ice  thickness.  Figure  4 1  shows  the  plots  of  v/fh  versus 
the  aspect  ratio (///;  and  the  velocity-to-thickness  ratio  vl 
h.  The  values  of  v/fh  are  in  the  range  of  0.18  and  0.24. 
This  means  that  the  average  movement  of  the  indentor 
during  each  ice  failure  varies  between  1 8  to  24%  of  the 
ice  thickness.  In  both  figures,  v/fh  remains  constant  as 
d/h  or  v/h  increase.  A  nondimensional  quantity,  fh/v 
(similar  to  the  Strouhal  number  for  fluids),  is  calculated 
and  listed  in  Table  3;  it  ranges  between  4.37  and  5.68  in 
this  test  series. 

Buckling  failure 

Tests  with  an  indentor  width  of  1 50  mm  and  an  ice 
thickness  30  mm  resulted  in  buckling  failure  (tests  N  N 
7 1  to  N  N  76  in  Table  3). 

Six  nondimensional  data  points  were  plotted  in  Figure 
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42,  in  which  theoretical  results  from  a  finite  element 
analysis  by  Sodhi  ( 1 979)  are  also  plotted  for  hinged  and 
frictionless  boundary  conditions.  The  ordinate  in  Fig¬ 
ure  42  is  a  nondimensional  buckling  pressure  FhlKL2(F 
is  ice  force,  d  is  the  structure  width,  K  is  the  specific 
weight  of  water  and  L  is  the  characteristic  length  of  the 
ice  sheet)  and  the  abscissa  is  the  ratio  of  the  indentor 
width  to  the  characteristic  length  of  the  ice  sheet.  The 
hinged  boundary  condition  means  that  there  is  no  rela¬ 
tive  displacement  between  the  ice  and  the  structure,  and 
the  frictionless  boundary  condition  means  that  there  is 
no  frictional  force  when  ice  moves  up  or  down  relative 
to  the  structure.  These  two  boundary  conditions  repre¬ 
sent  the  extreme  situations  possible  in  the  experiments, 
because  the  actual  boundary  condition  falls  between 
them. 

All  data  except  the  test  with  the  velocity  of  5  mm/s 
lie  in  the  zone  between  the  hinged  and  frictionless 
conditions.  With  decreasing  velocity,  it  seems  that  the 
nondimensional  buckling  pressure  decreases,  and  that 
the  boundary  condition  approaches  the  frictionless 
condition  from  the  hinged  condition. 
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Figure  41 .  Average  indentor  movement  per  ice  failure 
cycle  in  terms  of  ice  thickness  versus  d/h  and  v/h. 


d/L 


Figure  42.  Nondimensional  buckling  force 
F/dKL2  as  a  function  o/d/L  ( indentor  width 
d  =150  mm,  F  =  ice  force,  h  =  ice  thickness, 
K  =  9806  N/m  and  L  =  characteristic 
length ).  The  number  beside  each  data  point 
is  indentor  velocity  in  millimeters  per  sec¬ 
ond. 


SUMMARY 

We  conducted  92  indentation  tests  by  pushing  verti¬ 
cal,  flat  indentors  through  freshwater,  columnar,  float¬ 
ing  ice  sheets.  To  observe  the  failure  modes  and  to 
characterize  the  magnitude  and  nature  of  the  ice  forces, 
new  methodology  was  adopted  for  conducting  these 
tests.  For  instance,  forces  were  measured  by  supporting 
the  indentor  on  three  load  cells;  the  displacements  of  the 
carriage  and  the  indentor  were  measured  separately  to 
observe  the  transfer  of  energies  from  the  carriage  to  the 
structure  and  the  ice;  AEs  were  also  measured  and 
correlated  with  measured  ice  force.  The  indentor  width 
was  varied  from  50  to  1 50  mm,  the  relative  velocity  of 
the  indentor  from  1  to  9  mm/s,  and  the  ice  thickness 
from  20  to  60  mm. 

Instrumentation 

Although  indentation  tests  have  been  conducted 
earlier  by  other  researchers,  there  are  a  few  salient 
features  of  this  study  that  need  to  be  pointed  out. 

Ice  force  measurement 

The  interaction  ice  force  was  measured  in  two  ways: 

1 )  by  installing  a  load  cell  in  the  structural  support,  and 

2)  by  mounting  the  indentor  plate  on  three  load  cells. 
The  second  method  provided  the  actual  ice  force  that 
was  generated  between  the  ice  and  the  indentor  without 
being  influenced  by  the  deformation  and  vibration  of 
the  structure. 

Indentor  displacement 

Displacements  of  both  the  carriage  and  the  indentor 
were  measured  separately.  A  displacement  transducer 
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placed  on  the  floating  ice  sheet  measured  the  relative 
motion  of  the  indentor  with  respect  to  the  ice  sheet.  This 
allowed  us  to  compute  the  total  energy  supplied  by  the 
carriage  and  the  energies  dissipated  in  the  ice  and  stored 
in  the  structure. 

Acoustic  emissions 

Acoustic  emissions  were  recorded  during  the  tests, 
and  these  appear  to  be  correlated  to  the  microcracking 
activities  in  the  ice  and  the  ice  force.  The  cumulative  AE 
count  increases  with  increasing  ice  force. 

Test  results 

Microcracks 

At  the  onset  of  loading,  microcracks  immediately 
formed  in  front  of  the  indentor,  and  this  zone  of  micro¬ 
cracks  expanded  until  the  ice  failed.  During  high  veloc¬ 
ity  tests  ( v  >  4  mm/s),  we  observed  spurts  of  microcrack¬ 
ing  prior  to  each  ice  failure  event.  Microcracks  devel¬ 
oped  in  the  ice  ahead  of  the  indentor  during  the  time 
when  the  indentor  was  being  loaded,  and  then  the 
indentor  moved  forward  very  quickly,  crushing  or 
spalling  the  previously  damaged  (microcracked)  ice. 

Macrocracks 

Several  macrocracks  were  observed,  i.e.,  radial 
cracks,  circumferential  cracks  and  in-plane  cleavage 
cracks.  A  number  of  radial  and  circumferential  cracks 
nucleated  during  the  first  loading  and  failure  of  the  ice. 

Velocity  effects 

Depending  on  the  relative  velocity  of  the  indentor, 
ice  behavior  was  either  ductile  or  brittle.  Brittle  behav¬ 
ior  resulted  in  the  sudden  unloading  of  the  indentor, 
which  moved  back  to  its  original  undeformed  position, 
relative  to  the  carriage,  with  crushing  or  spalling  of  the 
damaged,  microcracked  ice.  Ductile  deformation  re¬ 
sulted  in  gradual  loading  and  unloading  of  the  indentor. 
It  appears  that,  even  fordifferent  indentor  velocities,  the 
indentor  displacements  at  the  time  when  the  peak  ice 
force  occurs  are  almost  the  same.  In  addition,  the  size  of 
the  damaged  ice  zone  during  the  first  loading  of  the  in¬ 
dentor  is  about  the  same  even  when  the  indentor  veloci¬ 
ties  are  different. 

Energy  analysis 

During  the  ice-structure  interaction,  energy  sup¬ 
plied  to  the  carriage  by  the  drive-screw  is  partly  stored 
in  the  structure  and  partly  dissipated  in  the  ice  as  a  result 
of  deformation  and  microcracking  in  the  ice.  During 
tests  at  high  indentor  velocity  (v  >  4  mm/s),  the  energy 
stored  in  the  structure  is  greater  than  that  dissipated  in 
the  ice.  This  means  that  the  deformation  of  the  structure 


is  larger  than  the  displacement  of  the  indentor  relative 
to  the  ice  sheet  until  the  ice  fails.  During  tests  at  low 
velocity  (v  <  3  mm/s),  the  energies  stored  in  the  struc¬ 
ture  and  dissipated  in  the  ice  increase  with  the  elapsed 
time,  but  the  energy  dissipated  in  the  ice  is  greater  than 
that  stored  in  the  structure.  This  means  that  thedisplace- 
ment  of  the  indentor  relative  to  the  ice  sheet  is  greater 
than  the  deformation  of  the  structure  during  the  first 
loading. 

For  tests  with  high  indentor  velocity  (r  >  4  mm/s). 
the  rate  of  energy  dissipated  in  the  ice  gradually  in¬ 
creases  with  the  elapsed  time,  and  the  curve  becomes 
steep  at  the  time  of  the  ice  failure.  For  tests  with  low 
indentor  velocity,  the  rate  of  energy  dissipated  in  the  ice 
is  much  greater  than  that  stored  in  the  structure,  and  it  is 
even  greater  than  that  supplied  to  the  carriage  after  a 
peak  in  the  force  record  occurs  and  the  structure  un¬ 
loads.  The  negative  values  for  the  rate  of  energy  stored 
in  the  structure  means  that  the  structure  gradually  swings 
back  to  its  original  position,  thereby  losing  its  stored 
energy. 

The  energy  used  to  produce  the  peak  force  is  almost 
the  same  even  fordifferent  indentor  velocities.  In  other 
words,  though  the  peak  force  and  the  displacement  of 
the  indentor  are  different  at  different  indentor  veloci¬ 
ties.  the  energy  dissipated  in  the  ice  up  to  the  first  peak 
in  force  is  almost  the  same  if  other  conditions  are  the 
same. 

Fortests  at  high  indentor  velocities  (v>4  mm/s),  the 
rate  of  energy  supplied  by  the  carriage  exceeds  that 
dissipated  in  the  ice;  the  rest  is  stored  in  the  structure.  As 
the  rate  of  energy  dissipated  in  the  ice  increases  with  the 
elapsed  time,  ice  is  weakened  by  microcracks,  and  the 
ability  of  the  ice  to  absorb  more  energy  decreases  with 
elapsed  time.  Just  before  the  peak  ice  force,  the  rate  of 
energy  supplied  to  the  ice  from  the  carriage  and  the 
structure  increases  rapidly  with  the  elapsed  time,  lead¬ 
ing  to  an  instability,  which  is  commonly  called  the 
brittle  failure  of  ice.  These  rates  are  much  lower  for  tests 
at  low  indentor  velocity  (v  <  3  mm/s)  than  for  tests  at 
high  velocity  and  attain  a  steady-state  value  as  opposed 
to  the  steep  rise  observed  in  high-velocity  tests.  The  low 
rate  of  energy  dissipation  in  the  ice  at  low  indentor 
velocity  leads  to  a  gradual  unloading  of  the  indentor  in 
a  stable  manner,  which  is  commonly  called  the  ductile 
failure  of  the  ice. 

Simultaneous  versus  lion-simultaneous  failure 

Because  the  pressure  transducers  were  damaged 
during  the  tests,  the  pressure  distribution  could  not  be 
determined  for  all  the  tests.  From  the  limited  data 
obtained  from  these  transducers,  we  found  that  a  nearly 
uniform  pressure  distribution  develops  during  the  ini¬ 
tial  contact  of  the  indentor  with  the  edge  of  an  ice  sheet. 
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In  a  few  tests,  the  interfacial  pressure  increased  to  a  high 
value  soon  after  the  indentor  contacted  the  ice  and  then 
remained  more  or  less  constant  while  the  ice  force 
increased  gradually  until  the  ice  failed.  More  tests  with 
wider  structures  need  to  be  done  to  determine  whether 
ice  fails  at  different  times  in  different  parts  of  the  contact 
area. 

Position  of  resultant  force 

In  many  of  the  tests,  the  interaction  forces  were 
measured  by  supporting  the  indentor  plate  on  three  load 
cells.  This  enabled  us  to  compute  the  position  of  the 
resultant  force.  We  were  looking  for  a  symmetrica! 
pressure  distribution  where  the  resultant  force  would 
remain  in  the  center  of  the  contact  area.  The  initial 
contact  point  between  the  indentor  and  the  ice  edge  may 
not  have  been  in  the  center  of  the  contact  area,  but  as  the 
load  increased,  the  position  of  the  resultant  force  moved 
to  the  center  of  the  contact  area,  indicating  a  uniform  or 
symmetrical  pressure  distribution  at  the  interface. 

First  peak  force  and  subsequent  peak  force 

Fi .  :»t  peak  forces  were  always  greater  than  the  subse¬ 
quent  peak  forces,  which  were  between  30  to  60%  of  the 
first  peak  force.  We  attribute  this  to  the  fact  that  ice 
offers  less  resistance  to  indentation  as  a  result  of  previ¬ 
ous  deformation  and  microcracking.  The  undamaged 
ice  offers  maximum  resistance,  thereby  resulting  in  a 
higher  interaction  force. 

Frequency  of  ice  force  failure 

The  frequency  of  ice  failure  increased  with  the 
increase  in  the  ratio  of  velocity  to  ice  thickness.  Using 
the  velocity  and  frequency  data,  we  calculated  the 
distance  of  indentor  movement  for  each  ice  failure;  it 
was  between  0. 1 8  and  0.24  times  the  ice  thickness. 

Buckling  failure 

The  ice  sheets  buckled  during  the  tests  when  the 
indentor  width  was  150  mm.  The  nondimensional 
buckling  load  was  found  to  be  between  those  predicted 
by  theoretical  results  for  frictionless  and  hinged  bound¬ 
ary  conditions.  Experimental  results  were  closer  to 
those  for  the  hinged  boundary  condition,  and  these 
decreased  with  decreasing  indentor  velocity. 

Scope  of  future  work 

The  data  obtained  in  this  study  provided  us  with 
many  interesting  results  that  were  only  possible  because 
of  the  new  methodology  used  to  measure  ice  forces  and 
the  measurements  of  carriage  and  indentor  displace¬ 
ment,  separately,  with  respect  to  a  fixed  datum.  There¬ 
fore,  we  recommend  that  all  future  experimental  studies 
of  this  nature  also  measure  forces  and  displacements  in 
a  similar  manner.  The  following  is  a  list  of  studies  that 


we  believe  should  be  conducted  for  a  better  understand¬ 
ing  of  ice-structure  interactions. 

1 .  Through  the  course  of  study,  our  observations  of 
the  deformation  of  the  structure  during  the  ice  loading 
showed  us  that  the  stiffness  of  the  structure  is  an 
important  parameter  during  ice  crushing.  The  frequency 
of  the  ice  failure  may  depend  on  the  structural  stiffness. 
Thus,  it  is  recommended  that  the  influence  of  the 
structural  stiffness  to  the  ice-structure  interaction  be 
investigated. 

2.  Observations  of  the  structural  response  at  the  time 
of  ice  crushing  are  needed  to  understand  the 
force-displacementrelationship. 

3.  Microcracking  and  energy/energy-rate  dissipa¬ 
tion  relationships  in  the  ice  should  be  studied  for  a  better 
understanding  of  the  ductile  and  brittle  behavior  of  the 
ice. 

4.  Further  studies  of  non-simultaneous  failure  in  the 
ice-structure  contact  area  should  be  conducted  to  deter¬ 
mine  the  influence  of  the  contact  area  on  the  magnitude 
of  the  ice  force. 
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Figure  45.  Test  NSC  <S'5. 

Figure  A6.  Test  NSC  S4. 
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Figure  A7.  Test  NSC  S3. 


Figure  AS.  Test  NSC  -S’. 
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Figure  A9.  Test  NSC  79. 

Figure  AID.  Test  NSC  78. 
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Figure  A 1 1  .Test  NSC  77. 


Figure  A12.  Test  NSC  76. 


Indentor  Width=100mm 
Ice  Th i ckness=5?mm 
Velocity”  ?,9inm/s 


Indentor  W ; dth= 1 00mm 
Ice  Th i ckness=5?mm 
Velocity”  ? . 9mm/s 


Time  ( s ) 

Time  ( s ) 

Figure  A13.  Test  NSC  69. 

Figure  A14.  Test  NSC  68. 
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Figure  A15.  Test  NSC  67. 


Figure  A16.  Test  NSC  66. 
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Figure  A17.  Test  NSC  65. 

Figure  A1S.  Test  NSC  64. 
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Figure  A 19.  Test  NSC  63. 


Figure  A20.  Test  NSC  62. 
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Figure  A21.  Test  NSC  61. 
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Figure  A22.  Test  NSC  60. 
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Figure  A23.  Test  NSC  59 

Figure  A24.  Test  NSC  49. 
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Figure  A29.  Test  NSC  44. 

Figure  A30.  Test  NSC  39. 
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Figure  A31.  Test  NSC  38. 
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Figure  A32.  Test  NSC  37. 
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Figure  A3  7.  Test  NSC  26. 

Figure  A3S.  Test  NSC  25. 
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Figure  A39.  Test  NSC  24. 


Figure  A40.  Test  NSC  23. 


Indentor  Width-100mm 
Ice  Th  i  ckness-57mm 
Velocity-  6 . 2mm/s 

i  _  . 


'F  /i 


Indentor  Width-100mm 
Ice  Th i ckness-5 1mm 
Velocity-  7 . Smm/s 


Indentor  Displacement  (mm) 

Indentor  Displacement  (mm) 

Figure  A41.  Test  NSC  22. 

Figure  A42.  Test  NSC  21. 
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Figure  A43.  Test  NSC  19. 


Figure  A44.  Test  NSC  18. 
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Figure  A47.  Test  N  N  31. 


Figure  A48.  Test  N  N  32. 
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Figure  A65.  Test  N  N  51 . 


Figure  A66.  Test  N  N  52. 
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Figure  A67.  TestNN  53. 


Figure  A68.  Test  N  N  55. 
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Figure  A69.  TestNN  56. 
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Figure  A71.  Test  N  N  58. 
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Figure  A70.  Test  N N 57. 
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Figure A72.  TestNN 59. 
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Figure  A73.  TestNN  60. 


Figure  A74.  Test  N  N  61 . 
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Figure  A75.  TestNN 62. 


Figure  A76.  Test  N  N  63. 
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Figure  A77.  Test  NN  64. 


Figure  A78.  Test  N  N  65. 
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Figure  A79.  Test  NN  66. 


Figure  A80.  Test  N  N  67. 
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Figure  A81  .Test  N  N  68. 
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Figure  A83.  Test  N  N  70. 
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Figure  A82.  Test  N  N  69. 
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Figure  A84.  TestNN  701 . 
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Figure  A89.  Test  N  N  75. 
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Figure  A90.  Test  N  N  76. 
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Figure  A91.  Test  N  N  81 . 
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supplied  to  the  carriage,  dissipated  in  the  ice  and  stored  in  the  structure  versus  time. 
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Figure  A92.  Test  N  N  31.  Figure  A93.  Test  N  N  33. 
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Figure  A94.  Test  N  N  35. 
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Figure  A95.  Test  N  N  37. 
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Figure  A96.  Test  N  N  39. 


Indentor  Wtdth»50mm 
Ice  Th 1 ckness»60mm 
Ve  1  oc t ty"6 . 9mmXs 


c-isai 


Indentor  Width-50mm 
Ice  Th  i ckness“60mm 
Ve  1  oc  1  ty”9 . 4 mm/'s 


r Indantor 


2 


••Structure 

/ICO 


Carr  t  ago 
loo 


2  3 

Time  ( s ) 


Figure  A97.  Test  N  N  41. 
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Figure  A98.  Test  N  N  43. 

Figure  A99.  Test  N  N  45. 
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Figure  A 104.  Test  NN  55. 
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Figure  A105.  Test  N  N  57. 
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